arXiv: 1507.0035Ivl [astro-ph.GA] 1 Jul 2015 


Mon. Not. R. Astron. Soc. 000, ^{5] (2015) Printed 3 July 2015 (MN style file v2.2) 


Fading Features Found in the Kinematics of the 
Far-Reaching Milky Way Stellar Halo 

Sarah A. Bird^’^ * and Chris Flynn^ 

^ Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, Chinese Academy of Sciences, 
80 Nandan Road, Shanghai, 200030, China 

^ Tuorla Observatory, Department of Physics and Astronomy, University of Turku, Vdisdldntie 20, FI-21500 Kaarina, Finland 
^ Centre for Astrophysics and Supercomputing, Swinburne University of Technology, Hawthorn, VIC, 3122, Australia 


Accepted 2015 June 26. Received 2015 June 25; in original form 2014 December 2 


ABSTRACT 

We test the long- term kinematical stability of a Galactic stellar halo model, due to 
iKafle et al.l (j2012[ l. who study the kinematics of approximately 5000 blue horizontal 
branch (BHB) stars in the Sloan Digital Sky Survey (SDSS). The velocity dispersion a 
and anisotropy parameter /3 of the stars have been determined as functions of Galac- 
tocentric radius, over the range 6 < i?GC < 25 kpc, and show a strong dip in the 
anisotropy profile at i?GC ~ 17 kpc. By directly integrating orbits of particles in a 
3-D model of the Galactic potential with these characteristics, we show that the a 
and P profiles quickly evolve on a time scale of a few x 10 Myr whereas the density p 
profile remains largely unaffected. We suggest that the feature is therefore transient. 
The origin of such features in the Galactic halo remains unclear. 
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1 INTRODUCTION 

Studying the Milky Way’s stellar halo is an important route 
to understanding galaxy formation, as the halo is such an 
old Galactic component. Intrinsically bright stars with eas¬ 
ily measured radial velocities have been the usual means of 
doing so, with red giants and horizontal branch stars as typ¬ 
ical tracers in such studies. Early studies of the stellar halo 
kinematics date to the 1950s, and focused on halo stars pass¬ 
ing through the Solar neighbourhood, but it was not until 
the 1980s that large (> 100) samples of halo stars tens of 
kpc from the Sun began to b e coll ected and an alysed (see 
the reviews bv ISanda^ (Il986l l and iHelniil ll2008ll l. 

Milky Way hal o BHB stars from ^ 5 to 50 kpc 
have b een studied bv ISommer-Larsen. Flynn fc ChristensenI 
(Il994h . They used about 100 stars to develop a kinemati¬ 
cal model of the outer Milky Way halo, with the surpris¬ 
ing result that the orbits of stars in the far outer halo 
(> 20 kpc) appear to be much more tangentia l than ra¬ 
dial. [Byjmj_SpmrneiMjaxsen_&^_Christens^ ([l9^) used sim¬ 
ulations of such stars orbiting in the Milky Way potential 
which showed such a distribution of halo orbits is stable over 
a Hubble time. 

Since then, nu merous studies have added to the sample 
of BHB halo stars lISommer-Larsen et al. 1 1 19971 : ISirko et al.l 


I 2 OO 4 I : iDeason. Belokurov fc Evcin3l201ll : iDeason et al.ll201^ 

but show a wide spread in the resulting kinematical models 
for the outer stellar halo. 

ISommer-Larsen et al.l (Il997l l analysed about 700 BHB 
stars, mainly within 20 kpc of the Sun, but also probing out 
to 50 kpc. They found that the outer stellar halo velocity 
dispersion (at « 50 kpc) was quite “cold” (i.e. low velocity 
dispersion), nearing 100 km s“^ compared with the value 
at the sun ~ 140 km s“^. They concluded that outer halo 
orbits must be quite tangential (with a tangential velocity 
dispersion of about 150 km s“^), given the observed density 
distribution of halo stars and assumptions about the Milky 
Way’s dark matter distribution. 

On the other hand, ISirko et al.l ll2004l l have advocated 
an isothermal outer halo {Rgc Pi Rq), in which all three 
components of the velocity dispersion are » 100 km s~^ , 
based on « 1200 BHB stars from SDSS. [Thom et ^ (l2005h 
subsequently analysed 530 BHB stars with radial velocities 
and distances from the Hamburg/ESO survey, finding it dif¬ 
ficult to discriminate between the simplest, isothermal kine¬ 
matic models and anything more complex, and advocating 
further studies of the inner halo to help resolve the issue. 


Very distant BHB stars have re cently been shown 
by Deason. Belokurov 8z EvansI (I2OIIII and IDeason et al.l 
(I2OI2I1 to have very “cold” kinematics - low velocity dis¬ 
persions of « 50 — 60 km s“^ in the radial range 100 to 
150 kpc. The density falloff in these regions is much steeper 
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than inside 100 kpc, and the dynamical times rather long, 
so the region is unlikely to be well mixed. Comparison with 
the kinematical models of the inner region, which implicitly 
assume the halo is well mixed, a r e thu s difficult. 

More recently, iKafle et al.l ll2012ll have used ~ 5000 
BHB halo stars found in the SDSS/Sloan Extension for 
Galactic Understanding and Exploration (SEGUE) to anal¬ 
yse the kinematics of the Galactic stellar halo. Work¬ 
ing from photometric distance estimates and radial veloc¬ 
ity measurements for each star, they performed a maxi¬ 
mum likelihood analysis to determine the velocity disper- 
sion a and anisotropy /3 profile, where jS is defined by 
iBinnev fc Tremaine! ll2008l l in spherical coordinates using ra¬ 
dial Or and tangential (ao, velocity dispersions such that 


/3 = 1- 


I 

2(jr 


( 1 ) 


iKafle et al.l (l2012f ) found a previously unseen feature, most 
prominent in the jS profile measured out to a Galactocentric 
radius of Rqc ~ 25 kpc, which shows a rapid decline at 
Rgc = 13 kpc, reaching a minimum at Roc = 17 kpc, 
followed by a sharp rise within just a few kpc (Fig. [T]). This 
feature has been confirmed in a very re cent study o f the h alo 
with an even larger sample of stars bv iKing et al.l (l2015l l. 

Given that the feature is so narrow and deep, but is 
made up of stars in an otherwise “hot” (i.e. high velocity 
dispersion) Galactic halo, we were motivated by the question 
of whether such a feature could be long-term stable or simply 
transient. 

We have tested this by setting up simulations of the 
stars in the Milky Way potential, and determined the 
model’s stability over the order of a Hubble time, Hnding 
that it is transient, and dissolves away in just a few tens of 
Myr. 

In Sec. [2] we d escribe our choice of potential to use 
for ou r simulations (iFlvnn. Sommer-Larsen fc Christensen! 
li^ )- updating it to be consistent with th e remarkable 
new c onstraints on the Galactic potential by IBovv &: Rbd 
(l2013tl . In Sec. [31 we describe our simulations testing stabil¬ 
ity over time of the stellar halo density distribution p, the 
velocity dispersion a, an d anisotropy P pr ofiles, motivated 
as such by the analysis of iKafle et al.l ll2012ll . We discuss our 
results and draw conclusions in Sec. (4] 


2 GALACTIC POTENTIAL MODEL 

We used a Milky Way potential mod e l sim ilar to that of 
iFlvnn. Sommer-Larsen fc ChristensenI (Il996l l. The poten¬ 
tial consisted of the sum of three components, namely the 
the potential due to the dark halo, a central component, 
and the disc. The dark halo potential was spherical with 
mass of order 10^^ Mq within 100 kpc. The central poten¬ 
tial was modeled by the sum of a spherical bulge/stellar- 
halo and a spherical inner core potential. The disc po- 
tential itself consisted of three Miyamoto-Nagai potentials 
llMivamoto fc Nagail 1197^1 . In this model, the disc has a 
scalelength of fl p = 2.2 kpc, to be consistent with the 
measurements of iBovv fc Rixl (l2013l l from the kinematics of 
over 16,000 G-type dwarfs in the SDSS/SEGUE survey dis¬ 
tributed between Galactocentric radii of 5 < i?Gc/kpc < 12. 
Our Galactic potential model parameters are listed in Ta- 
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Figure 1. Anisotropy (3 profile as a function of Ga lactocentric 
radius , Rqc-> showing /3 as calculated from the input [Kafle et alj 
ll2012h radi a l and tangential velocity dispersions, cfj. and at. 
iKafie et alJ l|2012l'l measured the anomalous feature at its ex¬ 
tremum to be = —1.2 at Rqq = 17 kpc. The error bars are 
Poissonian. 
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Table 1. Adopted Parameters for the Galactic Potential 


Component 

Parameter 

Value 

Unit 

Dark Halo 

rn 

15.0 

kpc 


M(Rgc < 100 kpc) 

~ 10i2 

Mo 


Uh 

220 

km 

Bulge/Stellar Halo 

rCi 

2.70 

kpc 


Me, 

3.0 

10®Mo 

Central Comp. 

rc2 

0.42 

kpc 


Me, 

16 

1O®M0 

Disc 

b 

0.3 

kpc 


r-Di 

5.81 

kpc 


Mn, 

106 

lO'^Mo 


r-Dj 

17.43 

kpc 


Md2 

-45.8 

IO^Mq 


r-Da 

34.86 

kpc 


Md3 

5.24 

lO'^Mo 


ble[T] The parameters have been set such that the rotation 
curve, i.e. circular velocity as a function of Galactocentric 
radius, is flat out to Rgc = 500 kpc — i.e. the outer limit 
of our modeled halo (although we only analyse the stars 
within 100 k pc). W e adopted the same parameters used by 
iKafle et al.l (120121 1. namely the Galactocentric position of 
the Sun at Rq = 8.5 kpc and the velocity of the local stan¬ 
dard of rest as ulsr = 220 km s“^. 


3 NUMERICAL METHODS AND 
SIMULATIONS 

We tested for stability over time of the Galactic stel¬ 
lar halo density distribution p, velocity dispersion cr, 
and anisotropy P profile s, setting them u p to initially 
match the observations of IKafle et al.l (l2012fl . We ran sim- 
ulations for 10 Gyr using 2.24 x 1 0^ pa rticles in the 
iFlynn, Sommer-Larsen fc ChristensenI (Il996l l Galactic po- 
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Figure 2. Radial density p profile of simulated halo stars. The 
initial dens i ty at 0 Gyr is the brok en power law a ssume d by 
[ 2012 !) (as obse rved by IWatkins et alj ||2009|) and 
Deason^ Belokurov Evans! ll201lh h where p oc Rqq and a = 2.4 
at Rqc ^ 27 kpc and ol = 4.5 at Rqc > 27 kpc (red crosses with 
red fit line). The profile is plotted at 0 and 10 Gyr (blue circles), 
showing the stability of p as we see negligible change between the 
initial and final states. The break radius at which the two power 
laws meet is Rqq ~ 27 kpc, or log(RGc) ~ 1-43 kpc. 
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Figure 3. U pper panel display s the initial velocity dispersion 
a profile from iKafle et all ll2012l ~). Radial and tangential velocity 
dispersions, cjr (blue crosses) and at (red asterisks) (km s“^), 
respectively, are plotted as a function of distance from the galactic 
centre, Rqc (kpc). The lower panel shows the a profiles after 2 
Gyr. The initial turnover features at Rgc ~ 17 kpc in the model 
halo disappear within 0.02 Gyr and the new a profiles remain 
stable in the simulations. 


tential as described in Sec. [2] The stellar velocities were 
initially randomly drawn from Gaussian distributions and 
set up in a spherical non-rotating configuration. The density 
profile (Fig . [2l) fo llowed the dou ble power law assume d by 


Kafle_et_^ (12^^, as observ e d bvl Watkins et al.l l|2nn9l) and 
Deason, Belokurov fc Evans! (1201 ll ) where o oc Rqc with 


a = 2.4 at Rqc ^ 27 kpc and a = 4.5 at Rqc > 27 kpc (red 
crosse s with red fit line). Very interestingly, iDeason et al.l 
ll2014h found recently a striking drop in the stellar density 
profile, with a = 6 at 50 kpc and a = 6 — 10 out to 100 
kpc. Although the stellar distribution of our simulations was 
truncated at 500 kpc and 0.1 kpc, we do not include the most 


recently found drop in stellar density as it would have neg¬ 
ligible effects on the stability of the system at much smaller 
radii near 17 kpc. 

The initial velocity dispersions profiles are shown in the 
upper panel of Fig. |3l 

Our simulations are of tracer stars in a fixed Galac¬ 
tic potential and the stars did not interact. The stellar or- 
bits were integrated using a Rung e-Kutta scheme control 
dPress, Flannery fc TeukolskvillOSB ) with a 0.01 Gyr (adap¬ 
tive) stepsize over a period of 10 Gyr, and the density p, 
velocity dispersion a, and anisotropy j3 profiles of the parti¬ 
cles as functions of Rqc were determined at 2 Gyr intervals. 

We found that the p profile was stable over the 10 
Gyr simulation, as seen in Fig. [5] when comparing the neg¬ 
ligible changes between the initial (red crosses) and final 
(bl ue circles) p r ofiles . Unlike the stable p profile, the ini¬ 
tial iKafle^^eTal] ll2012l) o profiles as seen in the top panel of 
Fig. |3]quickly relaxed to a new state, appearing in the lower 
panel of Fig. O which was stable over the remaining simu¬ 
lation time. Fig. ID shows that the initial profile was quickly 
changed within the first 0.02 Gyr (top right panel of Fig. [4]). 
By 0.1 Gyr (lower right panel of|4ll, relaxation had occurred 
and the new a profiles were stable over the remaining 8 Gyr. 
Similarly, the initial feature in /3 disappeared and the profile 
remained stable during the rest of the simulation. In addi¬ 
tion, we tested variations in the initial conditions which were 
consistent with th e boundaries set by the error bars of the 
iKafie et al.l ll2012l) determinations of sigma versus Rqc- In 
each test, the results remained similar, showing the outcome 
is not dependent on one unique choice of initial conditions. 

We experimented by initially exaggerating the feature 
to see if the stars would settle upon the actual model. We be¬ 
gan the simulation with a larger dip in the at profile, starting 
at a value of 50 km s“^ within the first 10 kpc instead of 90 
km s“^. All other conditions remained the same as our pre¬ 
vious simulations. We found that the cr profiles flatted and 
the anomolous feature dissolved away. We ran the simula¬ 
tions with the same conditions as previously described, but 
with a model potential with a perfectly fiat rotation curve , 
very similar to the one determined in iKafle et al.l ll2012l) . 
The results were the same as the other tests. 

Thus we concluded the Milky Way’s a and P profiles are 
in a transient phase, and the observed features at 17 kpc will 
smooth away within a few ten million years. 


4 SUMMARY AND CONCLUSIONS 


We tested the kinematic model of the Galactic halo proposed 


by K^e_e^_^ (|201^, by directly integ r ating particles in the 
iFlvnn, Sommer-Larsen fc ChristenserJ (Il996l ) 3-D model of 
the Galactic potential, under the assumption of an initial 
Gaussian velocity distribution. We found that the particles 
relaxed qui t e qui ckly from the initial profiles observed by 
iKafie et al.l (12012 ). The features at 17 kpc, most prominently 
seen in the anisotropy p profile as the steep decline reaching 
minimum at Rqc = 17 kpc and the sharp rise within a few 
kpc, dissolved after 0.02 Gyr. We propose that this feature 
in the Milky Way’s velocity dispersion cr and anisotropy P 
profiles is in a transient phase, and will flatten within a few 
ten million years. 

Observational studies have shown that stellar inhomo- 
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Figure 4. Similar to Fig. [S] but shown over the first 0.1 Gyr in the simulation. Radial and tangential velocity dispersions, dr (blue 
crosses) and at (red asterisks) (km s“^), are plotted as a function of distance from the Galactic centre, Rqc (kpc) . The <j pro fi les ar e 
plotted every 0.02 Gyr as indicated in the top left corner of each panel. The model input cr profiles, as observed bv iKafle et al.l ll2012t) . 
are plotted in each panel as the orange line (cr) and the green dashed line (ct). The dr and dt profiles change within the first 0.02 Gyr 
(top right panel). After 0.1 Gyr, the profiles have relaxed (lower right panel) and change negligibly over the remainder of the simulation. 


geneities and subst ructures exist in the halo, such as the 
Virgo overd ensitv d.Tur ic et al.l |2008| ) and Sagittarius stel¬ 
lar stream (llbata. Gilmore fc IrwinI 1 1994), which a r e pre - 
sumably remnants of merging satellites. Kafle et al.l (|2012| ) 
found that masking these two substructures did not signifi¬ 
cantly change their modeled proHles, although many under¬ 
lying stellar structures yet undetected are likely to exist, 
such that may or may not disrupt flat profiles and lead 
to temporary anomalous features. Such stellar structures 
are testa ble with simulations of halo formation around disc 
galaxies. IPeason et al.l (|2013h used results from such simu¬ 
lations to show that shell-like structures in the radial veloc¬ 
ities of halo stars are tightly related to the furthest orbital 
distance of stars left over from accreted satellites. 

Another possibility could be that the features 
in t he profiles r e flect underlying properties of the 
halo. ICarollo et al.l |2003) found a prograde inner halo 
and a retrograde more metal-poor outer halo; and 


ICarollo et al.l (120101) and Ide Jong et al.l ll2010l) found the 
transition between such halo components to lie between 
15 < fioc/kpc < 20, wh ich is the location of the feature 
in the iKafle et al.l ll2012h a and P proHles. In contrast to 
stellar substructures, th e presence of a dual halo has been 
deba te d for and aga in st (ISchonr ic h. Asplu nd fc Casagrandel 


201ll: Beers et al. [2O I2I: IPerinani fc Schonriclu 2013h 


Schonrich. Asnlund fc Casagrandel |2014^ . Despite the sug¬ 


gested explanations, the origins of (or combination there 
of) such transient features in the halo remain unclear. 

Alternatively, it is possible that the features in the cr 
and P proHles are due to the br oken power- l aw ste llar den¬ 
sity distribution p adopted by iKafle et al.l ll2012l) , as this 
function is marginalised over when determining the stellar 
velocity dispersions. We plan to repeat similar simulations 
to the current paper using t he more plau sible density of the 
alpha-beta-gamma profile of IZhad (Il996l ). which is a double 
power law but smooth. 
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In order to find what causes the features in t h e Gal actic 
kinematic profiles as interpreted by iKafle et al.l ll2012li , we 
plan to set up a complete library of dynamically stable halos 
in the adopted Milky Way potential. The library will be par- 
ticularlly useful to compare with the large amounts of new 
data anticipated from surveys such as Gaia and LAMOST. 
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